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If the Hartle-Hawking wave function is the correct boundary condition of our universe, the history 
of our universe will be well approximated by an instanton. Although this instanton should be 
classicalized at infinity, as long as we are observing a process of each history, we may detect a 
non-classicalized part of field combinations. When we apply it to a dark energy model, this non- 
classicalized part of fields can be well embedded to a quintessence and a phantom model, i.e., 
a quintom model. Because of the property of complexified instantons, the phantomness will be 
naturally free from a big rip singularity. This phantomness does not cause perturbative instabilities, 
as it is an effect emergent from the entire wave function. Our work may thus provide a theoretical 
basis for the quintom models, whose equation of state (EoS) can cross the cosmological constant 
boundary (CCB) phenomenologically. 
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I. INTRODUCTION 

One of the crucial tasks of quantum gravity is to understand the singularities in general relativity. When we 
consider the initial singularity of our universe [T] , the problem is related to various issues of physical cosmology, 
e.g., the origin of the emergence of time, the initial condition of inflation, the typicalness of our universe, etc. 
To deal with these issues, the traditional approach is to investigate the canonical quantization and to study 
the wave function of the universe [2]. 

After invoking the canonical quantization that includes the metric, what we eventually obtain is the master 
wave equation, the Wheeler-DeWitt equation. This equation is a partial differential equation and hence it 
requires boundary conditions. We do not know what should be the correct boundary condition, but perhaps 
the ground state of the universe can be a reasonable choice. Hartle and Hawking (HH) |3] suggested that 
the Euclidean path integral provides a good analog of the ground state wave function. For cosmological 
applications, the 0(4) symmetric metric ansatz would be a good simplification; and the Euclidean path integral 
can be approximated by the steepest-descent approximation, or by sum-over instantons. When we consider 
the Euclidean instantons, we need to complexify the time and hence every fields should be complexified by 
analyticity |1]. However, as long as the held is complex-valued, classical properties can never be restored in 
terms of equations of motion. Therefore, after the Wick rotation, the reality of the metric and the matter held 
is required: this is the classicality^ condition 


By ‘classicality’ we mean that a universe is classical, where the universe is originated from the wave function and the wave 
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Already some techniques have been investigated to calculate HH instantons and to estimate the probability 
distribution of each initial conditions mi- Typical expectations of the HH wave function are as follows: (1) 
it provides slow-roll inflation to obtain a classical history and (2) it does not prefer large number of e-foldings. 
The former is useful, but the latter is not a good news for inflationary cosmology [5]. However, if one considers 
more sophisticated models, this difficulty can be resolved. Note that when our universe begins, all field should 
be realized and satisfy the classicality condition. We envision that there exists some fields that in the early 
universe: an inflaton (or inflatons) that induces inflation, heavy mass fields and light mass fields compared with 
the inflaton. Regarding such a setting, the followings should be noticed. 

- If the mass scales of the fields are similar, then it is reasonable that these fields are equally excited at 
the same time. This may be related to the assisted inflation of multi-fields that can help to prefer large 
e-foldings [5]. 

- If there is a much heavier mass field (or fields), then in order to classicalize the heavier mass field, the 
lighter field should be excited uni. This excited lighter field can in principle be the inflaton field, which 
may further explain the preference of sufficiently large e-foldings. 

- Some modifications of the gravity sector in the early universe may help to prefer large e-foldings m- 

While this is not yet settled, it is fair to say that the HH wave function remains a reasonable theoretical basis 
for our inflationary universe |12j . 

If so, then the natural next question is, what will happen to the much lighter fields compared with the 
inflaton? Of course, at once they exist from the beginning, then these light fields should be regarded as a part 
of instanton. At the first glimpse, it is natural to assume that these light fields should be classicalized, too. 
However, if a field is decoupled from our phenomenological fields (standard model particles) and the amount of 
energy of this field is much smaller than that of the inflaton field, then even though the field is not classicalized, 
there is no way to distinguish the light field during and after the primordial inflation. As time goes on, however, 
the super slow-rolling and non-classicalized field can leave some distinguishable effects in the universe around 
the dark energy dominated era. This is a kind of ‘residue’ from the quantum gravity. Then can we see these 
effects in this universe? (Regarding this topic, for extended calculations, see [13] ■) 

Motivated by this scenario, in this paper we study the properties of a field that has negligible amount of 
energy compared to the inflaton, which is super slow-rolling and non-classicalized. By non-classicalized, we 
mean that the scalar field is not entirely realized from complex values (following the notion of classicality in 
13 El). Although this field has the form of a quintessence field, however, due to its non-classicalicity, some 
part of this field will also possess the phantom behavior. Therefore, effects of the non-classicalized field can 
be very well-embedded in a quintessence -I- phantom dark energy model, i.e., the quintom model [14] . This 


function itself is not classical. This is different from another notion of ‘classicality’ in the literature of inflationary physics; 
in this context, people consider a classicalization of quantum fluctuations. In this paper, our physical object and interest are 
different from the latter issue (classicalization of quantum fluctuations). 
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quintum model is known to be useful to investigate late time cosmology, especially in order to explain the 
crossing phenomenon of the dark energy equation of state over the cosmological constant boundary. Now the 
question is this: if there remain effects from a non-classicalized field as a quintum model, then what will be the 
signatures to our late time universe? This is the task of this paper. 

This paper is organized as follows. In SEC.|llj we briefly summarize previous results on the HH wave function. 
In SEC. |IIl| we discuss the behavior of the non-classicalized field that is indeed a quintom model; we also discuss 
the physical implications of this model. In SEC. EYl we discuss further interpretational issues, and finally, in 
SEC. 1^ we summarize this paper and discuss future issues that should be further investigated. 


II. HARTLE-HAWKING WAVE FUNCTION FOR TWO SCALAR FIELDS 


A. Basic formalism and classicality 


The ground state wave function by Hartle and Hawking [3] is defined as the Euclidean path integral for a 
compact 3-dimensional manifold S as a functional of the 3-metric and the field value x by 

o-Sb 


Jm 


( 1 ) 


where the 4-metric and the field (j) (for multi-field case, include all fields) take the value and y on 
E = dAi, where TW is a compact 4-dimensional Euclidean manifold. We integrate over all Ai that have E as 
their only boundary. 

In this paper, we investigate Einstein gravity with two minimally coupled scalar fields (we choose the 
units c = G = h = 


se = -J ^ . 


( 2 ) 


For the purpose of demonstrating qualitative properties, here we invoke a simple quadratic potential with mass 
mi and m 2 : 


( 3 ) 


V{(t)i,(t)2) = Vb -b ^mj(t>l + ^ml(l)l. 

What we want to attain are the following conditions: 

1. Vq is much smaller than m^: Vojmi <C 1. Therefore, during the inflationary era, we can ignore Vq- 

2. 02 satisfies over-damped conditions even with Vq'- rnd^jViii < dn (or, m2/H < 3/2, where = SttVo/S). 
Therefore, after the inflation era, 02 still satisfies the over-damped condition. 


a. Minisuperspace model We impose the minisuperspace model following the 0(4) symmetric metric ansatz 

( 4 ) 


j 2 + a^{T)dri^ 

uSj5 = ^ . 


mt 


From this choice of metric, it is convenient to redefine 


p = 


m 2 

mi 


( 5 ) 
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FIG. 1: Left: An instanton solution is defined on the complex plane t = x + iy. Right: By choosing a contour C 
(red arrows) we can draw a combination of the Euclidean and the Lorentzian manifolds. If we choose a proper initial 
condition and a proper turning time X, we can satisfy the classicality condition at large Y. 


The HH wave function is now 


where the action is reduced by (here, we ignored the Vq term) 

dr 


5e = ^ 


+ “) + 11 + 2 + '^ 2 ^ + '/’!+ 


( 6 ) 


(7) 


Even though /r ^ 1, we explicitly retain this term to study the behavior of the field (j> 2 - Along the contour C, 
the metric a starts from zero, which will be interpreted as the South Pole; along this contour, it grows to the 
boundary value b in the Lorentzian regime where (j)i takes the value Xi (FIG.[^. 

b. Steepest-descent approximation To approximately estimate the path-integral, we use the steepest- 
descent approximation. We approximate the wave function by summing over on-shell histories, the instan- 
tons, that satisfy the same boundary conditions [3]. For such an on-shell history p, the HH wave function is 
approximated by 


^[^Xi,X2] - X!® 

p 

Note that the on-shell condition is to satisfy the following equations of motion: 



where ' denotes a derivative with respect to r. 
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c. Classicality condition Since our universe follows the Lorentizian signature, a time contour in the path 
integral (Eq. should connect from Euclidean to Lorentzian manifold. The contour of r is defined on the 
complex plane (left of EIG.[^. The field values at the boundary of the scale factor b and scalar fields Xi should 
be real numbers. However, these metric and scalar fields are naturally complexihed along the complex time 
contour. We are interested in the condition of the endpoint (b and Xi)- By using the analyticity, we can choose 
a contour t = x + iy ior Q < x < X and 0 < y <Y (right of EIG. that connects from r = 0 to the endpoint. 
This contour connects from r = 0 to the turning point ai t = X through the Euclidean time; then, one can 
Wick-rotate to the Lorentzian time until the boundary at t = X + iY. 

If the action along a given history is complex-valued and if the real part and the imaginary part of the action 
rapidly vary up to the variation of canonical variables, then the Hamilton-Jacobi equation is not satished and 
hence the history is no more classical. On the other hand, if the real part of the Euclidean action varies slowly 
compared to the imaginary part, then the Hamilton-Jacobi equation (the classical equation of motion) will be 
approximately satisfied. According to [6], this is called the classicality condition: 


|V^ Re S'e[^,Xi,X 2 ]| < |Va Im Ssib,Xi,X 2 ]\, 


( 12 ) 


where A = b,Xi- practice, the classicality condition can be presented by 

llm a| \lm (bj\ 

^ (El^ ^ ^ 

for all f’s as t increases, and hence correspond to the reality at the endpoint m- 

When the classicality condition is satisfied, we can interpret that the instanton generates a universe along 
the time direction. For a classical universe, one can approximate the probability of the Wheeler-DeWitt wave 
function by 

P[b,Xi,X2] a \nb,xuX2f - (14) 


d. Initial conditions The boundary condition at the South Pole comes from the regularity condition, 

a{T = 0 ) = 0 , a' {t = 0 ) = 1 , = 0 ) = 0 . 

At the end endpoint, we impose the following conditions where b and Xi s-re real values: 

a(r = X + iY) = b, (f>i{T = X + iY) = Xi- 


(15) 


(16) 


At the turning time, because of the analyticity, we impose the Gauchy-Riemann condition: 

= = = ?^ir = X)Ar = X). 

OX toy ox lay 


(17) 


This system is constructed by second order differential equations of three complex-valued functions: a and 
(j>i. We have eight boundary conditions at r = 0 and three conditions at the endpoint. We solve this problem 
by choosing a scalar field value at r = 0, 


= 0) = (/>z(0) = |0i(O)|e*^' 


(18) 


where |0i(O)| and 9i are real. One can solve this initial value problem to calculate time evolutions of a and (j)i 
from r = 0. For a given |(()i(0)|, in order to satisfy classicality conditions, one needs to tune X and 9i. 
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B. Summary of previous results and motivations 


e. Applications to primordial inflation These conclusions are already proven by previous authors: 

- For a single held inhaton with V = Vq + {l/2)mf(fA , if mf IVq < Gtt and hence if the potential is in the 
slow-roll regime, then the probability distribution is consistent with that of the quantum held theory in 
de Sitter space m- 

- On the other hand, if rnf > Gtt, then (j) cannot be classicalized around 0 = 0. This was proven 
analytically as well as numerically in UU- 


-Asa simple extension, if there are two helds 0i and 02 with V = (l/2)TO^0f + ( 1 / 2 )to 20 | and mi/m 2 S> 1 
(hence, 02 direction is a slow-rolling direction), then to classicalize the heavy mass direction 0 i around 
01 = 0 , we must require the condition uni 


■ ^ „ ■ j - < Gtt. (19) 

(l/ 2 )mi 0 i 

This in turn requires 02 ^ (mi/m 2 ) to classicalize both helds (and, this initial condition is the most 
probable one as well, see details in [ID]). 


In the early universe, there may exist various helds. To classicalize heavy helds, some slow-rolling helds need 
to be excited and these excited slow-rolling helds can be the origin of inhation. 

/. Motivations: what about a slower direction? If the inhaton held is excited, inhation is turned on, and 
as the inhaton decays, matters and structures will be formed. However, what will happen if there was a much 
slower direction than the inhaton held? Let us call this held a quintessence. 

If this quintessence is decoupled from the other matter helds and its direction rolls much more slowly than 
the inhation itself, then even though the held is not classicalized, it would not induce any observable effect. 
Hence, even though the quintessence field is not entirely classicalized, during and post inflation, it renders no 
observable impact. 

However, at late times after radiation and matter dominant eras, such quintessence held may in principle 
exhibit some physical imprints. Then what will be the signatures of the non-classicalized quintessence to our 
late time universe? This physics should be connected to physics of the dark energy, which is the task of this 
paper. 


III. PHYSICS OF NON-CLASSICALIZED FIELD: QUINTESSENCE AND/OR PHANTOM 

We explicitly write the relevant quantities as a = Or- + zoi, 0i = (fir + icfu, and 02 = 02 r + * 02 i- Let us 
assume that mi ^ m 2 , where 0i is the infiaton field and 02 is the quintessence field. In addition, let us further 
assume that a and 0i are almost completely classicalized, while 02 is not. In other words, as t —> Y (where 
t = Y is almost the end point of inflation), 

' ' ’ 101.1 ’ 102 . 1 “ 


a, 


0 ( 1 ). 


( 20 ) 





In addition, we further assume that around the turning time t = X, 


I « 1, (21) 

so that the contribution from (j )2 does not affect inflationary physics (hence, when the universe is created, the 
probability is mainly determined by (j)i and does not sensitively depend on ^ 2 )- 

Then after the primordial inflation and matter/radiation dominated era, there will be an era dominated by 
the quintessence held. 



A. Behavior of non-classicalized field 


g. Equations of motion Equations of motion in Lorentzian signatures are separated by real parts 
0 = dr + ^'^1- ^'^Ir + ^'ir ~ <Pli ~ ~ 2 (‘(’ir ~ *(*1* + ^2r ~ 4'2i) 

-Oi ^20ir<^li + ‘2^2r4>2i ~ 4>lr4’li ~ 4’2r4’2^ i 


Stt 

0 = ‘(j)\r 3 

0 = ^2r 3 

and imaginary parts 


Qj'pQj'p ~\~ CL'iCL'i 

al + af 
CLrpGj>p I 
a'i + af 


1 

-3| 

( —OrOi + diar \ 

V Or + af ) 

) 02r 

-3| 

( —drUi + CliGr \ 

V a2 + a| ) 


4>li + 4>lr, 
4>2i + 


0 = 


8 7 ^ / \ 

di + + ‘^4>2r4>2i ~ 4‘lr‘fli ~ 'p2r4>2ij 

+ ^<^lr + ^2r ~ 4>\i — 4>2i~ 2 ~ ^2r ~ 4' 


0 — 4^ii H“ 3 
0 


Qi'fCL'p ~\~ CL'iOj'i 

al + aj 


. Qjj'di -\- didr . 

Li + 3 1 - — -^- I (pir + 01i, 


. drdr H” didi . 


a; + a‘ 

d'fd'i ~\~ d'ld'p 


4>2r + M ^2z 


( 22 ) 

(23) 

(24) 


(25) 

(26) 
(27) 


where ' is with respect to the Lorentzian time. 

h. Existence of ai,(j)ii —>■ 0 turning time We hrst argue that there exists a turning time when Ui —>■ 0 and 
4>ii —>-0. As a toy model, let us hx = 0.01, |(/'i(0)| = |^2(0)| = 0.9. Since ^ 1 and the initial held 
position is the same, the total energy contribution is dominated by (fi. 

If we classicalize two helds at the same time, then the optimized point is 9i ~ —0.1676, 02 — —0.0016; and 
along the turning time X ~ 0.85, we obtain the classicalized Lorentzian history. We can solve the same initial 
condition not only along one time contour, but also over the complex plane (FIG.[^ (see also [IS]). This result 
shows that along the turning time X ~ 0.85, three curves (dashed, dotted, and thin white curves, corresponding 
Oi = 0, (fii = 0, and (j) 2 i = 0, respectively) coincide and hence along the Lorentzian time, all helds will be 
classicalized. 

Now let us consider the situation that we tilt 02 from the optimized value and violates the classicality of 4)2 ■ 
As long as 1 / I ^ 1) the effects of (j )2 will be very restricted for a and (fi. If the tilted angle increases, 
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Re a Im a 



FIG. 2: flr, ai, 4)\r, 4 ’' 2 -r, and (j) 2 i over the complex time plane t = x + iy, for = 0.01, |0i(O)| = |<;^2(0)| = 0.9 

with initial conditions 6 \ ~ —0.1676, and 62 — —0.0016. Dashed, dotted, and thin white curves are ai = 0, cj>ii = 0, and 
4 > 2 i = 0 , respectively, where we superimposed three curves on the figure of middle-right. 


then by tuning a proper 6 * 1 , again we can obtain a good turning time where —?► 0 and (/)ii 0 are satisfied. 

For example, in FIG. we tilt 62 and check that there still exists a turning time X that satisfies a^, —>■ 0. 

For more realistic applications, in FIG. we demonstrated a case when the tilted value is much larger 
02 = 9o — 24007r/8192 to demonstrate a phantom phase. In this case, we choose y = Q to apply for a realistic 


cosmological model that should satisfy (see details in SEG. IIIB). Even though the tilted value is larger 
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FIG. 3: (f>ii for slightly tilted 02 = ± 27r/8192, where 6o is the optimized value ~ —0.0016, the plus sign is for left, 

and the minus sign is for right. Dashed, dotted, and thin white curves are at = 0, (f>ii = 0, and (f) 2 i = 0, respectively. 


than the optimized value, still the classicality of a and (j)i is robust. 

i. Embedding in quintom model FIG. has shown the existence of a history that satisfies > 0 and 
(pii —>■ 0. We have already demonstrated this numerically. To be prudent, we further check its consistency 
through analytic approximations. In this regard, if we choose the proper turning time that approximately^ 


satisfies —>■ 0, 4>n —>■ 0, and hr jar = H, then equations are simplified by 

0 = H + H'^ + — - 'p2i~ 2 '^2r ~ 

0 = ^ir + (^9) 

0 = 2p2r^2i — fJ-'^4’2r4’2i, (30) 

0 = (f>2r “f 'iH(j)2r + (31) 

0 = 02i + 3iF02i + E^'p2i- (32) 


Therefore, except Eq. (30) that is related to a^, this system of equations are indistinguishable to the system of 
a quintessence field and a phantom field. 


We already found that there exists a direction that satisfies —>■ 0 and hence Eq. (30) should be consistent 
in the end. We can further check the consistency. During the inflationary regime, we can approximate H as a 
slowly varying function. Then the follows are solutions: 


(j)2r = A+e “+‘ + A_e (j)2i = B+e “+* + i 3 _e 


(33) 


where 


a± 


3H ± y/9H^ - AgB 
2 


(34) 


^ In numerical analysis, |ai|/|ar| and |i?iii|/|0ir| rapidly approaches to zero and hence (although and (pu are not exactly zero) 
this is a very good approximation. 
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FIG. 4: Ur, ai, (jiir, <^ii, 4 ’' 2 -r, and 02i over the complex time plane t = x + iy, for = 0, |</)i(0)| = |02(O)| = 0.9 with 
initial conditions 6 \ ~ —0.1676, and 62 = 0o + A 2 , where A 2 = —24007r/8192. Dashed and dotted curves are a; = 0 and 

(i>ii = 0 . 


If we insert this to Eq. (301, then this term behaves as 


(35) 


and hence as time goes on Eq. (301 will be satisfied. This implies that as time goes on, i.e., as ai and (jin decay 
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02;m/02r[t] 



FIG. 5: 02i/(^2r for 62 = Oo — 24007r/8192 after the turning time X, where 9o is the optimized value ~ —0.0016. Here, 
each color denotes different (black: 0.01, red: 0.001, blue: 0.0001). The ratio approaches a constant as time goes on. 


to zero, Eq. (301 will be automatically satisfied^. 


j. Initial conditions From the above analysis, we thus have various possible initial conditions in the post¬ 
inflation period. Let us discuss them in the following: 


If A 4 . — B- = 0, then 


10 : 


2i\ 


102 




0 , 


(36) 


and hence the classicality of 02 is satisfied. In other words, the classicality of 02 is only allowed by a 
finely-tuned initial condition. 


- If A± and B± are all non-zero, then 


I02d 

I02r| 


B_ 

—> — = const. 


- If A_ = Bj^ = 0, then 

102.1 

and hence the phantom dominance. 


(37) 


(38) 


If there exists an instanton from the natural parameter space, then |02i|/|02r| const is the most reasonable 
condition. Of course, for realistic calculations, H varies with time and hence details are quite complicated. 
However, as long as we consider the time when a and 0i are sufficiently classicalized, still this assumption 
I02i|/|02r| const is quite reasonable from numerical calculations. In FIG.[^ we show that the ratio |02i|/l02r| 
approaches to a constant as time goes on for various choices of 

In conclusion, this model is well embedded in a model with a quintessence field ipq and a phantom field ijjp 
with the initial conditions satisfying IV'pl/IV'gl — const. 


® We may further choose H = /r/ 2 to automatically cancel Eq. j30| , but we will not further restrict our parameters. Since 
\o,i/ar \ <K 1, as we see in Eq. \22\ , contributions to energy-momentum tensors will be well approximated by a quintessence field 
and a phantom field. 
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B. Generalization: implications to late time cosmology 


k. Generalization of Hartle-Hawking inspired quintom model According to the above analysis, for a scalar 
field system 


S' D y d'^x^/^ 


(</>!, </> 2 ) 


and after the classicalization of metric and inflaton field (fi, at the end of inflation, it can be transcribed to a 
two-field model as: 


S D y d*Xy/^ 




where 


— 02r; 

Ipp — 4^2i 5 

U{tpq,tpp) = Re = 0, 02r +# 2 i)- 


(39) 


(40) 

(41) 

(42) 


This model, with the signs of the kinetic energy of the two fields being opposite from each other, is actually the 
“quintom” model [14]. One salient property of this model is that its EoS can have a crossing behavior around 
the cosmological constant boundary w = — 1, either from above to below or vice versa. There exist varieties of 
ways to realize quintom behavior (for earliest realizations, see [T7| for two-field models and [TS] for single field 
models with higher derivative). In this paper, we realize a quintom model in a more fundamental way, i.e., 
from the Hartle-Hawking wave function. 

1. Implications to late time cosmology In our example, we choose the potential to be of quadratic form: 


U{ipq,'ijjp) = m|Re[(((> 2 r + Z(/>2i)^]/2 = m'^i'ipq — t/’p)/2. Then according to action (39l, the total energy density 
and pressure of this quintom model are 


P=\i’q- - V^p), P=^'>Pl- ^V’p - - V'p), 


and the equations of motion for ifq and ifp are: 


(43) 


(44) 


(45) 


ifq + SHijjq + m\ipq = 0, 'ifp + SHifp + m^'tpp = 0, 

respectively. Thus the equation of state of the whole system is: 

_ - V'p - ml - ifl) + 2pi + 2po 

i’l-i’p+ ml {ipl - ipl) + 2pi -H 2po ’ 
where pi and pi are contributions from (jji] pq and po are contributions from the cosmological constant Vq- 
Let us focus on the following points, which has been shown in FIG. [^ 

- During the inflationary phase, the contribution of Vq was negligible. However, after inflation ends, pi and 
Pi become negligible, while Vq may eventually emerge. 
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y=oo 


-H 

y - 

w —> -1 

We are observing here 


w < -1 


Vo begins to be dominated. 

i 

Primordial inflation finished. 


(|)i classicalized, |(|) 2 i|/|(|) 2 r|~ 0 (l) 

,1 ► » 

"- 


FIG. 6: A diagram of the conceptual history. 


- In this limit, (j )2 remains in the over-damped regime, since we assumed to^IVq < Gtt. Then 

4>2r,2i OC e 5 

where H is determined by Vq. We see, once again, that 

IV'pI 

r = ~ const. 

m\ 

This ratio r will be determined when the field is created by an instanton. 


( 46 ) 


( 47 ) 


In this limit, from Eq. (45) we have 


%-%- - V'p) - Wo 


1 -I- w = 


- i’l) 


V’g-V’p + wi(V'2-V'p) + 2Fo’ ' - il^l) + 2Vo 


( 48 ) 


If we set the initial conditions such that < '4’p, then it is natural to have 1 -|- tc < 0, i.e., the phantom 
behavior. However, along with the evolution, the field energy density will eventually become negligible relative 
to the constant term Vo, and the EoS will approach the cosmological constant boundary (CCB) w = —1. To 
see this, it is useful to define the energy density and pressure for each field component as: 

^ ]2 , ^ 2 i2 1/2 I 2/2 1;2 l2;2 1/2,12/2 /-^inN 

Pq= 2 Vq +Pq = Pp = -Pp = -l^Vp + i^ru^Vp. ( 49 ) 


such that Pq > 0, Wq = Pq/Pq > —I, Pp < 0, Wp = Pp/Pp > —I. Furthermore, from the equations of motion one 
gets Pq ~ ar IppI ~ 3(1-1-™,)^ both pq and the absolute value of pp decrease with time. This means 

that both il>q and tpp will have decreasing contribution in the universe, while Vq remains a constant. This is 
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FIG. 7: Left: Evolution of the equation of state w with respect to Ina^ in our model, where is the classicalized 
scale factor of our universe. Right: Evolution of the energy density fraction of dark energy Q.de with respect to In a,, 
in our model. In the numerical study, we choose m 2 = 8.0 x 10“®^Mp (black), 5.0 x 10“®^Mp (red), 1.0 x 10“®^Mp 
(blue) respectively, while Vo = 0.5 x 10“^^®Mp. Initial conditions: ipqi ~ 0.33Mp (black), 0.58Mp (red), 3.13Mp 
(blue), ijjpi ~ —0.005Mp (black), 0.055Mp (red), 0.68Mp (blue), 'ipqi — 3.67 x 10“®^Mp (black), 2.79 x 10“®^Mp (red), 
1.55 X 10“®^Mp (blue), ippi ~ 4.93 x 10“®^Mp (black), 4.80 x 10“®^Mp (red), 4.56 x 10“®^Mp (blue). 


why the universe will eventually be dominated by Vq, having w approaching —1. However, since the evolution 
of the two fields are the same except for the initial condition, the relation between ipq and tpp could be more 
subtle. If during the evolution it happens that ipp exceeds '0^, then w will become larger than —1, and the 
quintom behavior will appear. 

In FIG. we draw three cases of evolutions in our model. We start from a phantom phase with w < — 1, 
with different initial conditions. One can see from the plot that although the initial values are different, they 
all eventually converge to the w = — 1 line, which confirms the above analysis. Moreover, two of the three lines 
display crossing behavior, and the other one approaches —1 directly from below. We also plot the evolution 
of the energy density fraction VIoe for the three cases. All of which shows that in the future VIee —>■ — 1, 
namely the universe will be dominated by dark energy. Actually, all the other components (including ipg, tpp, 
matter, radiation, etc) decays other than the constant term Vq, so it is an attractor solution that the universe 
will always be dominated by Vq. Furthermore, our plot shows that at the current time (Ina = 0) we have 
w ~ 1.1, rioE — 0.68, which are well within the newest Planck data, which suggests that w = —(2cr, 
Planck2015 TT+lowP)^ and LIa = 0.686 ± 0.020 (la, Planck2013) [H]. 

One important remark is that this model can also be free from the big rip singularity. According to pp] . 


^ From joint analysis of data, the best fitted value of w could be closer to —1, for example w = —l.OO60Q Qg® based on Planck 
power spectra, Planck lensing, and external data m- 
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when the universe is dominated by the dark energy with w, the time scale At for the universe to be of size a 
is approximately 


At = 


3(l+u;)ffo%/l-n„o 
1 m „ 


/ 3(1 + ^) \ 


UI ^ - 1 , 
w = —1, 


(50) 


Ho — 

where Hq is the current Hubble parameter and flmo is the current density fraction of matter in our universe. 
A big rip singularity occurs when a —>■ oo, which will cause: 


{ At —)■ oo ti; = — 1, 

(51) 

^ ~ 3{l+w)HoVl-n^o W < -1. 

Since in our scenario when dark energy dominates the universe {floE 1), w already always converges to 
(or larger than) —1, and therefore it must correspond to the condition that At goes to infinity. That is, it is 
impossible for the big rip singularity to occur in a finite time in the future. 

In summary, through an explicit example, we showed that our Hartle-Hawking instanton solution can be 
applied to late time cosmology, with the light fields behaving as phantom and quintessence fields in the quintom 
model. Since there exist a future attractor where w = —1, the big rip singularity is also avoided. Thus one 
may say that Hartle-Hawking interpretation of the quantum universe can also provide a theoretical basis for 
the quintom dark energy models, whose EoS can cross the CCB. 


IV. INTERPRETATIONS 

The ground state wave function can be represented by the Euclidean path integral [3] 

^'o [ht,u,x]= J e"'®"". 

This Euclidean analytic continuation is the origin to introduce complexified fields. The necessity to introduce 
complexified fields is very clear from some examples, by comparing calculations using instantons and using 
quantum field theory in de Sitter space m- These complexified fields are not a problem in general, since we 
require the reality at the endpoint of the path integral (e.g., asymptotic future infinity). 

However, a problem appears in our universe, since we are not at the endpoint but in the process. If we are not 
seeing the exact endpoint^ then it is allowed to see some effects of the imaginary part of a field, i.e., a ghost-like 
behavior of a scalar field. Since the instanton approximates this wave function, it already contains quantum 
contributions. Hence, the instanton and its imaginary part are an emergent result of the entire path integral. 

Can we find an analog of this phenomenon? Hawking radiation can be an example. Hawking radiation can 
be interpreted by using a particle propagator [21]. The particle propagator can be approximated by a classical 
path over the Euclidean analytic continuation. This process can be interpreted as follows: a particle comes 
out from the event horizon backward in time (or oppositely, one can say that a negative energy particle comes 
into the black hole forward in time) and the same energy particle is detected at the asymptotic future infinity. 
The classicality will be imposed at the future infinity; but as long as a particle moves backward in time, the 
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bulk description cannot be classical. Now if we cut a Cauchy surface including inside the event horizon, the 
Cauchy surface includes ghost-like particles. This can be conceptually related to the fact that the renormalized 
energy-momentum tensor {T^u) can violate the null energy condition around the horizon. Even though the null 
energy condition is violated, it does not cause a serious instability, since the effects of the negative energy are 
emergent results from the entire path integral. 

Of course, there are some conceptual differences between black hole physics and cosmology. For a black hole 
case, the renormalized energy-momentum tensor is an averaged result (7)jjy), not an independent instanton. 
On the other hand, for a cosmological case, we are in a special universe and hence we should see a special and 
independent instanton. Can we justify this phenomenon further? We remain this for a future work. However, 
in conclusion, it seems that if our universe could be phantom-like (i.e., w < —1), this Hartle-Hawking inspired 
quintom model can be a legal way to justify phantomness in terms of quantum physics. 


V. CONCLUSION 

In this paper, we investigated the Hartle-Hawking wave function with a two-scalar-field model. This wave 
function is well approximated by summing over instantons. In general, these instantons will be complexified, 
but in order to obtain a well-defined probability, one needs to require the classicality of each instanton, i.e., all 
fields should be realized at infinity. However, as long as we are an observer not at infinity but at a finite time, 
it is permissible to observe the imprints of the imaginary part of the fields. 

In order to embed this possibility to the late time cosmology, we assumed two massive canonical scalar fields 
{4>i is an inflaton and (j )2 has a slower direction) plus a cosmological constant with some physical conditions 
imposed: (1) initially the energy contribution of is dominant over (j )2 and (2) after decays, (j )2 still satisfies 
the over-damped condition. Then during primordial inflation, the scale factor a and the inflaton field (j>i will 
be realized sufficiently; and as long as the first condition is satisfied, even if (j )2 is not realized, the realization 
of a and (j>i can still be robust. 

Then all effects of (j )2 can be negligible during inflation; as our universe approaches the dark energy dominated 
era, however, the non-classical and super-slow-roll scalar field will contribute to the equation of state. If the 
amplitude of the imaginary part of (j )2 is larger than that of the real part of then w < — 1 can be attained. 
However, as time goes on, all real and imaginary parts must decay to zero, and hence the EoS will either cross 
the cosmological constant boundary lu = — 1 then reduce to it, having a quintom-like behavior, or go to — 1 
directly like phantom models. In either of the two ways, the EoS only stays below —1 for a finite time, so there 
should be no concern about the big rip singularity problem. Thus our model has shown that Hartle-Hawking 
wave function can be viewed as a theoretical basis and a possible origin of the quintom dark energy models in 
late time cosmology. 

Usually, the phantomness can be easily introduced by a ghost field. However, a ghost field causes perturbative 
instability, and hence physically disallowed [22] . In this paper, the imaginary part of a scalar field behaves as a 
ghost field with negative kinetic energy; but this term came from a non-perturbative effect of the entire wave 
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function. Therefore, we may say that this phantomness can be an emergent effect of quantum gravity. 

In this paper, we only restricted to quadratic potential, but in principle it can be generalized to various 
potentials based on different motivations. In addition, one may apply the same philosophy to investigate other 
physical phenomena such as black holes. If further investigations can indeed establish the connection between 
dark energy and the non-classicallized instantons, then this would be the first evidence of effects emergent from 
quantum gravity. 
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